Hamiltonian renormalization for bound state problem in gluodynamics

The systematic approach to study bound states in gluodynamics is presented. The method utilizes flow equations together with low-energy phenomenology, that provides the perturbative renormalization scaling in conjuction with the change of the basis to constituent gluon states. The renormalized to the second order effective Hamiltonian of gluodynamics is obtained at low energies, which provides the constituent description for glueball bound states. The approach allows to include perturbative QCD corrections into nonperturbative calculations of manybody techniques. The performed numerical calculations support the constituent picture of hadron observables.